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摘 要 : 明晰 土壤 盐分 空 
羊 河 下 游民 勤 盆 地 为 研究 区 
局 的 分 析 ,探讨 民 勤 盆地 土壤 盐分 空 
(1) 实测 土壤 盐分 含 
势 且 值 域 范 围 随 士 层 加 次 而 变 大 , 正 ( 负 ) 空 

反 )。(3) 变异 函数 均 以 指数 模型 为 最 好 , 尼 随 十 
化 规律 一 致 , 随 着 土 层 加 深 空 


,通过 对 0~30 cm 土 层 
间 异 质 性 


量 和 变化 范围 均 随 土 层 加 深 而 变 小 
间 自 相关 怕 


区 人 研究 
ARID ZONE RESEARCH 


间 变 异性 是 土壤 盐 渍 化 程度 研究 和 盐碱地 改 
对 土壤 盐分 原 位 数据 统计 学 特征 A 
变化 规律 ,以 其 
, 且 均 呈 偏 正 态 分 布 。 


E 均 随 距 离 的 增 大 而 减 小 (0~10 cm 负 的 空 
层 加 深 而 变 大 (0.62 一 0.69 一 0.81) , RSS FIFE 0.0002 , CJ (Co+C) Fl AoE 
s 间 异 质 性 由 弱 到 强烈 再 到 中 等 (4) 插值 精度 10~20 cm > 20~30 cm > 0~10 cm。3 
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本 研究 选取 干旱 内 陆 河流 域 石 
空间 自 相 关 性 .变异 性 和 分 布 格 


良 的 前 提 。 


为 研究 区 土壤 盐 渍 化 治理 提供 理论 依据 。 结 果 显 示 : 


(2) Moran”s 了 指数 均 星 ~” 型 波动 趋 
s 间 自 相 关 性 相 


个 土 层 土壤 盐分 均 存在 空间 分 异 特性 
(不 同 土 层 ) 上 ,分布 规 律 由 简单 到 复杂 


非 盐 化 土 (<1.30%) 和 4 种 不 同 程度 盐 渍 土 在 3 个 土 层 中 均 


20.00% ) > 轻 度 ( 约 9.00%) > $R E( < 0.50%). 
关键 词 : 土壤 盐分 ; 空间 变异 ; 地 统计 学 ; 


土壤 盐 涡 化 是 自然 环境 和 人 类 活动 共同 作用 


,水 平方 向 (同一 土 
,从 聚集 状态 (0~10 ecm) 逐渐 趋 于 斑 块 化 (10~20 cm) 和 斑点 化 (20~30 cm). 


层 


) 上 , 盐 渍 化 程度 均 自 西南 向 东北 逐渐 增加 ;垂直 方向 


9 分布 , 盐 渍 土 面 积 中 度 ( >70.00% ) > 强度 ( 约 


民 勤 盆地 ; 石 羊 河 下 游 


导致 土壤 盐分 累积 或 散失 量 的 差异 “"。 通 过 对 比 不 


影响 下 的 土壤 退化 现象 ,也 是 全 球 性 的 生态 环境 问 
题 之 一 ,不 笃 济 的 发 展 ,也 危害 了 生态 
环境 的 稳定 性 "。 土 壤 盐 分 一 般 在 表土 层 (0~30 cm) 
分 布 较为 广泛 ,其 量化 后 可 以 表征 盐 渍 化 的 进程 ， 
同时 盐分 在 不 同 土 层 方向 的 分 布 与 移动 导致 了 盐 
渍 化 的 纷繁 多 变 ””。 而 在 盐分 运 移 过 程 中 受 母 质 、 
环境 和 生物 等 因素 影响 而 具有 空间 异 质 性 ,并 通过 
空间 上 斑 块 之 间 casa guinea ai 进而 呈现 出 
不 同 斑 块 之 间 的 ( 自 ) 相 关 性 5”。 通 过 获取 空间 异 质 
性 与 其 自 相 关 性 有 利于 土壤 资源 的 科学 利用 和 明确 
盐 渍 化 斑 块 对 环境 的 影响 。 因 此 ,通过 量化 土壤 盐 


同 土 层 之 间 的 盐分 变化 可 以 掌握 不 同 土 层 中 盐分 的 
差异 ,有 助 于 更 清楚 地 了 解 土壤 盐 渍 化 进程 ”。 综 上 
所 述 ,结合 土壤 盐分 量化 和 土 层 之 间 的 空间 分 异 特 
征 可 以 更 深入 地 理解 盐 渍 化 现象 。 

民 勤 盆地 作为 干旱 区 荒漠 -绿洲 生态 系统 的 典 
型 代表 ,风沙 . 盐 清和 干旱 使 境内 土壤 盐 活化 较 严 
重 ”。 现 有 研究 多 从 地 下 水 、 地 表 水 的 水 体 矿 化 度 

出 发 中 ,到 土壤 盐分 离子 的 变化 ，, 以 及 不 同 植被 


群落 的 二 盐分 与 养分 的 空间 变化 “” ,而 缺乏 对 
民 勤 盆地 土壤 盐分 在 空间 上 的 系统 分 异 研究 。 本 


文 以 土壤 盐分 量化 为 基础 , 辅 以 表土 层 不 同 层次 之 


分 来 获取 区 域 化 的 盐分 数据 ,可 以 在 不 同 区 域 中 展 
示人 研究 区 内 盐 渍 化 的 空间 分 布 ,并 对 土壤 盐 演化 动 
态 监 测 具 有 重要 意义 。 同 时 , 细 化 表土 层 后 ,研究 发 
现 盐 分 受到 结构 因素 和 随机 因素 的 综合 影响 ,从 而 
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1 材料 与 方法 


1.1 研究 区 概况 

民 勤 盆地 (38?24'~39?15'N, 102°40'~103°55’ E) 
位 于 甘肃 省 武威 市 民 勤 县 ,处 于 石 羊 河流 域 下 游 且 
由 腾 格 里 沙洲 和 巴 丹 吉林 沙漠 三 面 ( 除 西 面 ) 包 围 
(图 1)。 属 温带 大 陆 性 干旱 气候 区 ,四 季 分 明 ,干燥 
少雨 多 风 , 年 均 温 7~8 % ,风沙 频 发 ,年 均 风 速 2.7 
ms ', 降 水 量 110 mm 左右 ,年 平均 获 发 量 2483 
mm。 作 为 流域 的 尾 间 段 ,沙漠 侵袭 和 成 水 灌溉 导致 
土壤 沙化 、 盐 化 加 剧 。 土 壤 类 型 主要 是 灰 棕 漠 填 、 
沙土 和 盐 化 土 。 植 被 以 典型 的 沙 生 、 旱 生 和 盐 生 植 
物 为 主 , 白 刺 (Nitraria spp.) H48 (Haloxylon ammo- 
dendron) 和 树 柳 (Tamarix chinensis ) 是 其 主要 的 优势 
灌木 。 
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图 1 研究 区 地 理 位 置 及 采样 点 分 布 


Fig. 1 Location and distribution of sampling points 


in the study area 


12 土壤 盐分 数据 采集 及 测定 

选取 人 为 因素 影响 较 小 且 植 被 分 布 相对 均匀 
的 区 域 ,每 个 采样 点 以 南北 方向 间隔 100 m 布设 10 
mx10 m 样 方 ,3 个 记 为 重复 。 样 方 内 南北 对 角 线 取 
3 个 点 的 土 样 (深度 为 0~10 cm、10~20 cm 和 20~30 
cm) 并 分 层 混 匀 ,四 分 法 装 袋 后 带 回 实验 室 备 用 。 
2020 年 (79 个 样 点 ) 和 2021 年 (51 个 样 点 ) 共 计 130 
个 采样 点 ,经 自然 风干 REAR OP L a Bos 
准 土 样 用 于 盐分 测定 。 利 用 样品 电导 率 及 与 烘 干 
残渣 法 联合 获取 盐分 含量 ""。 参 考 《 中 国 盐 涡 土 》 
和 《土壤 农 化 分 析 》 盐 涡 化 类 型 划分 指标 对 人 研究 区 
土壤 盐 渍 化 进行 分 类 分 级 (<2 ge kg SEER AG ; 2~ 


3 geko FARES EE 3~5 oko 中 度 盐 渍 十 ;5~10 gkg 
强度 盐 渍 土 ;>10 g'kg 盐 澳 土 ) ” ,并 评价 盐 渍 化 
程度 。 

1.3 研究 方法 

1.3.1 正 态 性 检验 及 方差 分 析 Kolmogorov-Smirnov 
检验 (K-S 检 验 ) 是 一 种 常用 的 拟 合 优 度 检验 方 
法 ,可 根据 样本 数据 推断 其 来 自 的 总 体 是 否 服从 
某 一 特定 理论 分 布 ,在 5% 置 信 水 平 下 进行 土壤 盐分 
数据 的 正 态 分 布 检验 ,并 根据 P 值 大 小 得 出 结论 "1。 
单 因素 方差 分 析 一 般 用 于 两 个 及 以 上 样本 平均 数 
差别 的 显著 性 检验 ,用 于 研究 某 一 个 类 型 的 自 变量 
与 男 一 个 数值 型 的 因 变量 之 间 的 差异 (5%)"。 

13.2 空间 自 相 关 性 FH Moran’ s 1 指数 展现 土 
培 盐 分 的 空间 依赖 性 ( 自 相 关 性 ) ,计算 公式 为 : 


Moran ’s [= 
"$ $w, -os / SSW Se O 
i=1j= i=1j= i= 


式 中 :n 为 需要 估计 的 样 点 数目 ;x 和 % 分 别 为 变量 x 
在 样 点 ;和 7 处 的 观测 值 ; x 为 x 的 平均 值 ;W; 为 相 邻 
权重 。 值 域 为 [-1, 1], 当 三 0 时 ,表示 变量 在 此 空间 
分 布 上 无 相关 性 , 且 呈 随机 分 布 ; 当 / 值 趋 于 1 时 , 变 
量 之 间 的 空间 分 布 关 系 呈 现 完全 正 相 关 和 聚集 状 
态 。 反 之 ,7 值 趋 于 -1 时 ,变量 呈 完 全 负 相 关 且 分 散 
分 布 。 
1.3.3 空间 异 质 性 利用 半 变 异 函 数 的 拟 合 参数 解 
析 土 壤 盐分 分 布 的 空间 异 质 性 特征 ,并 从 随机 性 和 
结构 性 变异 程度 两 个 方面 进行 解释 。 选 取 块 金 值 
(Nugget, Co) 4 8A (Sill, C+C)、 偏 基 台 值 (Partial 
Sil,C) 和 变 程 (Range,4o)4 个 主要 参数 。Co+C 表 示 
变量 的 总 变异 性 ,其 中 Co 表示 随机 因素 引起 的 空间 
异 质 性 ,C 表 示 结 构 ( 系 统 ) 因 素 引 起 的 ;二 者 之 比 
Cu(CotC), 即 基底 效应 展现 由 空间 自 相 关 的 部 分 引 
发 的 空间 变异 程度 ;4 表示 变量 具有 空间 相关 性 的 
阅 值 大 小 , 即 在 阔 值 范围 内 该 变量 具有 空间 相关 
性 ,超出 范围 则 独立 无 相关 性 *”。 
1.3.4 Kriging 插值 ”Kriging 插 值 是 地 统计 学 中 评估 
空间 分 布 的 一 种 方法 ,包含 多 种 算法 ,常用 的 有 普 
i Kriging、 简 单 Kriging 和 泛 Kriging o KWAK 
用 普通 Kriging 法 对 土壤 盐分 进行 插值 。 人 公式” 为 : 
S50= 5,xW, (2) 


式 中 ;5 为 待 估算 盐分 含量 的 样 点 ;5; 为 已 知 样 点 i 
的 盐分 含量 ;到 为 样 点 ;的 权重 值 ;* 为 参与 计算 的 
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已 知 盐 分 含量 样 点 的 数量 。 
1.3.5 模型 精度 检验 ”为 了 保证 预测 值 的 无 偏 性 ， 
采用 平均 误差 (Mean Error, ME) .平均 标准 误差 
(Mean Standardized Errors, MSE) 、 标 准 均 方 根 预 测 
误差 (Root- Mean- Square Standardized Errors, RMS- 
SE) 、 均 方 根 预测 误差 (Root-Mean-Square Prediction 
Error, RMSE) 平均 标准 误差 (Average Standard Er- 
ror,ASE) 和 决定 系数 ( 尼 )6 个 指标 对 Kriging 捕 值 结 
果 进 行 精度 检验 。 
1.4 数据 分 析 

运用 GS+9.0 进 行 半 变 异 函 数 模 型 第 选 和 Mo- 
ran’ s 1 指数 分 析 , Kriging 插 值 以 及 空间 分 布 图 使 用 
ArcGIS 10.3 建 模 和 绘制 ,利用 Excel 2016 进行 统计 
分 析 。 


2 结果 与 分 析 


2.1 统计 分 析 及 正 态 性 检验 

人 研究 区 内 土壤 表层 十 上 层 (0~10 em) 中层 
(10~20 cm) 和 下 层 (20~30 ecm) 的 盐分 含量 呈 逐 渐 递 
减 的 趋势 , 即 随 士 层 变 深 盐 分 含量 均值 逐渐 减少 
(5.44 g.kg 一 5.23 g*kg'—5.07 g- kg), HERZL 
范围 缩小 (1.19~32.08 g- kg 一 1.24~19.05 g+ kg 一 


= 


25 p (a) 0~10 cm 
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1.43~16.29 g.kg )。 同 时 ,3 个 土 层 中 土壤 盐分 含量 
的 中 位 数 均 小 于 其 均值 , 且 标 准 差 为 3.50+0.38 g- 
kg ,说 明 50% 以 上 的 样 点 盐分 含量 在 均值 水 平 以 下 
且 离 散 程度 较 大 。K-S 正 态 性 检验 (95% 置 信 区 间 ) 
结果 表明 :0~10 cm 10~20 cm 和 20~30 cm 3 个 土 层 
土壤 盐分 经 对 数 变换 后 的 P 值 分 别 为 0.2902、 
0.7952 和 0.5170, 均 > 0.05, =Œ Mm Æ (Mm FE 0.2076~ 
0.8057) 正 态 分 布 (图 2)。 方 差分 析 结 果 表 明 ,不 同 
土 层 盐分 无 显著 差异 (P=0.164) ,因此 可 用 地 统计 学 
法 探索 其 空间 分 异 特征 。 
2.2 空间 自 相 关 性 

3 个 土 层 土壤 盐分 的 Moran's7 指数 均 呈 “人 一" 型 
波动 趋势 (图 3), 即 随 着 距离 的 增加 ,Moran s7 指数 
在 0 值 上 下 波动 , 且 正 空间 自 相 关 性 大 于 负 空 间 自 
相关 性 ; 值 域 分 别 为 [-0.151, 0.302] 、[ -0.185， 
0.392] 和 [-0.259, 0.552] ,范围 随 土屋 加 深 而 变 大 ， 
空间 自 相 关 性 和 空间 差异 性 均 随 土 层 加 深 而 变 大 ; 
除 0~10 cm 负 的 空间 自 相 关 性 ( 随 距离 增 大 而 增 大 ) 
外 ,其 他 土 层 的 正 ( 负 ) 空 间 自 相关 性 均 随 距离 的 增 
大 而 减 小 。 
2.3 空间 异 质 性 

通过 对 线性 球状、 指数 和 高 斯 4 种 变异 函数 模 
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图 2 样 点 盐分 含量 正 态 分 布 


Fig. 2 Normal distribution of salt content in samples 
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3 土壤 各 层 盐分 Moran 's 7 指数 


Fig. 3 Moran’s 7 index of soil salinity in each layer 


型 模拟 的 土壤 盐分 结果 进行 反复 比较 ,筛选 出 最 优 
半 变 异 函 数 模型 及 其 插值 相关 参数 ( 表 1)。3 个 土 
层 均 以 指数 模型 为 最 优 ( 尼 随 土 层 加 深 而 变 大 
(0.62 一 "0.69 一 0.81) ) , H. RSS (5&2 ) F492: 0.0002, Hi 
明 3 个 土 层 的 模拟 结果 均 可 以 有 效 表征 各 土 层 土 壤 
盐分 的 空间 异 质 性 。Co( 块 金 值 ) 在 0.008~0.060, 是 C 
FY 13.80% ~66.67% ,说明 在 模拟 过 程 中 随机 因素 得 
到 了 很 好 的 控制 ,较为 稳定 而 不 影响 模拟 结果 。Cw 
(Cot) (基底 效应 ) 随 土 层 加 深 由 弱 到 强烈 再 到 中 
等 ,说 明 盐 分 含量 在 不 同 土 层 中 的 空间 异 质 性 受 随 
机 性 和 结构 性 因素 的 影响 不 同 。4( 变 程 ) 的 变化 与 


因素 的 影响 相对 较 弱 ,以 致 盐 分 含量 变化 较 其 余 两 
层 剧 烈 ,空间 异 质 性 最 弱 量 4 最大。 

2.4 空间 分 布 格局 

2.4.1 插值 及 精度 验证 ”对 插值 结果 (图 4) 进 行 精 
度 验 证 ( 表 2) 可 知 ,3 个 土屋 ME( 平 均 误 差 )(-0.004~ 
0.006) 与 MSE( 平 均 标 准 误差 )(-0.018~0.017) 均 趋 
近 于 0, 筛 选 的 变异 模型 满足 对 土壤 盐分 预测 的 准 
确 性 与 有 效 性 ;ASE( 平 均 标 准 误差 ) 与 RMSE( 均 方 
根 预 测 误差 ) 差 值 的 绝对 值 均 <0.03, 对 变异 程度 的 
预测 控制 较 稳 定 可 信 ;RMSSE( 标 准 均 方 根 预 测 误 
差 ) 均 趋 近 于 1 ,对 变异 程度 的 预测 在 合理 范围 内 。 


Cu(Co+C) 一 致 。 其 中 ,0~10 cm EER ZNE 、 风 
沙 践踏、 晨 殖 等 随机 性 因素 的 影响 较 大 ,而 结构 性 


表 1 土壤 各 层 盐 


的 变异 


因此 ,参数 设置 合理 ,插值 精度 均 符合 要 求 , H 10~ 
20 cm > 20~30 cm > 0~10 cm, 


函数 理论 模型 及 相关 参数 


Tab. 1 Variogram theoretical model and related parameters of soil salinity in each layer 


ne KEEN 块 金 值 基 台 值 基底 效应 变 程 决定 系数 残 差 
Lelen 变异 模型 Co CC Col(Cot€) Aoki R RSS 
0~10 指数 模型 0.060 0.076 0.79 44.64 0.62 0.0002 
10~20 指数 模型 0.008 0.066 0.13 8.79 0.69 0.0002 
20~30 指数 模型 0.026 0.065 0.40 18.09 0.81 0.0002 
(a) 0~10 cm (b) 10~20 cm (c) 20~30 cm 


Bi | o 
非 盐 化 轻 度 中 度 强度 盐 十 


图 4 不 同 土 层 的 土壤 盐分 空间 分 布 


Fig. 4 Spatial distribution of soil salinity in different soil layers 
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表 2 不 同 土 层 土壤 盐分 插值 模型 验证 结果 
Tab.2 Verification results of interpolation model of soil salinity in different layers 

P 平均 误差 均 方 根 误差 标准 均值 误差 标准 均 方 根 误差 平均 标准 误差 
go ME RMSE MSE RMSSE ASE 

0~10 —0.004 0.25 -0.018 1.13 0.22 

10~20 0.001 0.22 0.002 0.98 0.23 

20~30 0.006 0.21 0.017 1.03 0.21 


24.2 空间 分 布 特征 3 个 土 层 土壤 盐分 在 空间 分 
布 上 存在 一 定 差异 性 (图 4), 非 盐 化 土 和 4 种 不 同 程 
度 盐 渍 十 在 3 个 土 层 中 均 有 分 布 。 水 平方 向 上 , 盐 
渍 化 程度 均 自 西南 向 东北 逐渐 增加 ,与 石 羊 河流 向 
一 致 ;垂直 方向 的 分 布 规律 由 简单 到 复杂 ,从 聚集 
到 分 散 。 研 究 区 98% 以 上 区 域 为 不 同 程度 的 盐 渍 
土 ( 表 3) ,以 中 度 盐 渍 土 为 主 ,面积 均 > 70.00% H. 
20~30 cm > 0~10 cm > 10~20 cm; 强度 盐 渍 士 次 之 ， 
主要 分 布 在 腾 格 里 沙漠 边缘 (研究 区 东北 - 北 - 南 沿 
线 边 缘 ) 区 域 , 面 积 随 士 层 加 深 而 减 小 (由 0~10 cm 
的 23.43% 减 少 到 20~30 cm 的 15.49% ) 且 逐渐 趋 于 
斑 块 化 ; 轻 度 盐 渍 土 位 列 第 三 ,主要 分 布 于 研究 区 


三 角 洲 为 例 ,通过 研究 土壤 盐分 空间 变异 的 变 程 指 
标 确定 出 每 1000 km? tie KEE 107 个 样 点 2 ,而 不 同 
地 点 有 着 不 同 的 样 点 分 布 数量 标准 ;从 采样 点 属性 
来 看 ,不 同 的 间距 与 协 变 量 (生境 、 地 表 和 覆盖 、 利 用 
方式 等 ) 的 选择 对 区 域 盐分 含量 估算 有 影响 全 ;从 
处 理 方式 来 看 ,干燥 与 否 .过 得 大 小 .测定 方法 选取 
等 不 同 的 处 理 都 会 对 盐分 含量 拟 合 模型 精度 有 影 
H; 从 结果 获取 来 看 , 拟 合 (数学 归纳 、 统 计 建 模 
等 ) 与 反 演 (卫星 .雷达 无 人 机 以 及 影像 的 光谱 选 
取 ) 都 会 给 结果 带 来 误差 的 。 因 此 ,在 样品 选取 与 
采集 过 程 中 需要 综合 考量 方法 . 协 变量 处理、 结 
展现 等 的 影响 ,才能 使 得 研究 更 加 准确 并 符合 
实际 。 


西南 毗邻 武威 例 地 区 域 , 随 着 土 层 加 深 研 究 区 内 部 
也 逐渐 出 现 斑 块 化 分 布 区 域 ,面积 变化 规律 与 强度 
Eb Yet LAA AR IR (3.12% 12.20% 14.15% ) ; $k E 
面积 最 小 (< 0.40%). TP ERE EE EEE 
西南 部 的 盆地 边缘 ,20~30 em 土 层 占 1.27% ,面积 最 
大 ;0~10 em 土 层 次 之 ,为 20~30 cm 土 层 的 12;10~ 
20 cm 十 层 最 小 , 仅 占 0.03%。 


3 讨论 
3.1 采样 点 质量 与 盐分 研究 精度 的 关系 


盐 土 的 盐分 组 成 与 母 岩 的 类 型 和 成 分 有 密切 
的 联系 ,为 了 获取 更 具有 代表 性 .典型 性 的 样品 , 根 
据 研 究 区 的 盐分 状况 .植被 类 型 等 确定 优质 的 采样 
点 2 。 采 样 点 的 确定 需要 结合 实际 ,综合 运用 多 种 
评价 准则 来 分 析 获 取 。 从 采样 点 数量 来 看 ,以 黄河 


3.2 沙化 对 盐 渍 化 程度 空间 异 质 性 的 影响 

不 同 土 层 不 同 斑 块 表现 出 空间 异 质 性 ”。 盐 
省 化 程度 分 布 基本 呈现 灌区 向 四 周 加 深 的 趋势 , 恰 
与 两 大 沙漠 侵袭 民 勤 盆地 的 方向 相反 ,土壤 沙化 程 
度 可 能 对 盐 渍 化 程度 有 影响 ”。 插 值 结果 表明 , 靠 
近 腾 格 里 沙漠 区 域 的 盐 省 化 程度 高 于 靠近 巴 丹 吉 
林 沙 漠 区 域 (图 4) ,实测 原 位 数据 结果 亦 是 如 此 ( 表 
4) ,原因 在 于 巴 丹 吉林 沙漠 的 侵袭 程度 要 高 于 腾 格 
里 沙漠 二 ” ,在 基于 出 人 沙 地 、 人 员 配置 等 因素 考量 
设置 腾 格 里 沙漠 边缘 区 域 的 采样 点 较 多 日 分 布 广 
泛 , 而 布设 在 巴 丹 吉林 沙漠 边缘 区 域 的 采样 点 相对 
较 少 旦 集中 分 布 在 沙漠 边缘 两 端 ( 图 5), 也 会 影响 
两 大 沙漠 边缘 区 域 土壤 盐分 的 空间 异 质 性 。 其 潜 
在 分 布 与 实际 情况 之 间 复 杂 的 变化 与 运 移 方式 有 
待 进一步 研究 ”。 


表 3 不 同 土 层 的 土壤 盐分 面积 占 比 


Tab.3 Proportion of soil salt area in different soil layers at different times 


土 层 /em 非 盐 十 /9% 轻 度 盐 溃 十 /9p 
0~10 0.63 3.12 
10~20 0.03 12.20 
20~30 1.27 14.15 


中 度 盐 涡 土 /% 强度 盐 渍 土 /% 盐 渍 土 /% 
72.44 23.43 0.37 
71.59 22.59 0.38 


75.80 15.49 0.09 


1020 干旱 区 研究 40% 
表 4 沙漠 边缘 区 域 统 计 信息 对 比 
Tab. 4 Comparison of statistical information of desert edge area 
区 域 采样 点 数 范围 Mg'kg) ”平均 值 /(g*kg") 标准 差 /g'kg) 标准 误差 Mg.kg) ”中 位 数 /(g kg’) 
巴 丹 吉林 沙漠 边缘 区 22 个 1.31~32.08 3.96 2.05 0.44 3.33 
民 勤 贫 地 腹地 70 个 1.29~23.18 4.86 3.77 0.49 3.82 
腾 格 里 沙漠 边缘 区 38 个 1.18~8.89 7.16 6.85 1.11 5.70 


图 5 沙漠 边缘 采样 区 域 的 对 比 


Fig. 5 Comparison of sampling areas at desert edge 


3.3 民 勤 岔 地 盐 渍 化 程度 的 变化 趋势 探究 

民 勤 盆地 是 基于 石 羊 河流 域 的 农业 垦殖 区 ?” 。 
20 世 纪 70 年 代 以 前 ,自然 条 件 下 的 蒸发 积 盐 是 该 地 
区 土壤 盐 渍 化 的 主要 原因 ”;70 年 代 之 后 ,大 量 地 
下 水 开采 灌溉 引起 的 土壤 盐分 人 为 再 分 配 加 剧 了 该 
地 区 的 次 生 盐 溃 化 “。 陈 丽 娟 等 和 顾 梦 知 等 ”的 
研究 发 现 , 青 土 湖 地 区 在 石 羊 河流 域 综合 治理 生态 


分 标准 (间隔 不 同 ) 也 会 对 土壤 盐分 空间 异 质 性 的 
表达 造成 影响 ,综合 多 个 土 层 的 土壤 剖面 研究 ”与 
仅 针 对 表土 层 甚 至 其 他 单一 土 层 的 结果 不 同 且 存 
在 一 定 差 距 ”“ 。 误 差 与 样本 量 对 于 试验 结果 的 影 
响 是 造成 空间 异 质 性 差异 的 又 一 影响 因素 ,包括 方 
法 选择 `、 人 为 误差 ,数据 累积 和 迭代 程度 以 及 系统 
变化 等 ”。 
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通过 筛选 最 优 变异 模型 获取 研究 区 土壤 盐分 
的 空间 分 布 格局 ,表现 不 同 土 层 盐分 的 空间 差异 ， 
结论 如 下 : 

(1) 3 个 土 层 土壤 盐分 含量 和 变化 范围 均 随 土 
层 加 深 而 变 小 , 且 均 呈 偏 正 态 分 布 ;Moran’'s 1 指数 
均 呈 “型 波动 趋势 , 值 域 范围 随 士 层 加 深 而 变 大 ， 
空间 自 相关 性 均 随 距离 的 增 大 而 减 小 。 

(2) 土壤 盐分 的 变异 因数 均 以 指数 模型 最 佳 ， 
模拟 结果 均 可 有 效 表征 各 土 层 土壤 盐分 的 空间 异 
质 性 有 旦 随机 性 因素 控制 较为 稳定 ;空间 异 质 性 随 土 


Ee 


恢复 工程 中 ,初期 (2010 年 ) 与 中 期 (2016 年 ) 以 强 
度 盐 渍 土 为 主 “ 汪 ,而 本 研究 发 现 , 工 程 结束 后 
(2020—2021 年 ) 中 度 盐 渍 已 达 研 究 区 的 73.55% ,说 
明 近 15 a 来 ,研究 区 盐 渍 化 程度 有 所 改善 ,由 强度 逐 
渐 向 中 度 转换 ,间接 表明 了 石 羊 河流 域 综合 治理 生 
态 恢 复工 程 对 民 勤 盆地 土壤 盐 渍 化 的 恢复 也 具有 
一 定 效果 。 
3.4 空间 异 质 性 变化 的 探究 

诸多 学 者 对 表层 土壤 的 盐 涡 化 程度 变化 提出 
了 不 同 的 思考 。 不 同 环境 下 ,生物 与 非 生物 因素 的 互 
作 关 系 是 影响 盐分 分 布 或 聚集 多 朝 的 关键 因素 ””。 
民 勤 地 区 因 所 处 环境 限制 ,植被 所 延伸 的 丰富 度 ” 
和 生物 量 所 等 因素 导致 了 不 同 土 层 变化 程度 的 差 
异性 ;地 下 水 和 和 土壤 湿度 所 等 水 分 条 件 又 会 从 非 
生物 的 角度 抑制 或 促进 这 些 变化 。 不 同 的 土 层 划 


层 加 深 由 弱 到 强烈 再 到 中 等 ,0~10 cm 土 层 因 受 雨 
雪 、 风 沙 . 战 踏 、 尾 殖 等 随机 性 因素 的 影响 较 大 , 空 
间 异 质 性 最 弱 且 变 程 最 大 。 

(3) 经 交叉 验证 插值 精度 (10~20 cm > 20~30 
cm > 0~10 cm) 均 符合 要 求 。 插 值 结果 表明 ,3 个 土 
层 土壤 盐分 均 存 在 空间 分 异 ,水 平方 向 (同一 土 层 ) 
上 , 盐 渍 化 程度 均 自 西南 向 东北 逐渐 增加 ,与 石 羊 
河 走向 一 致 ;垂直 方向 (不 同 土 层 ) 上 ,分布 规律 由 
简单 到 复杂 ,从 聚集 状态 (0~10 em) 逐渐 趋 于 斑 块 
化 (10~20 cm) 和 斑点 化 (20~30 em) 

(4) 非 盐 化 土 和 4 种 不 同 程度 盐 涡 土 在 0~30 
cm 土 层 中 均 有 分 布 , 以 中 度 盐 渍 土 为 主 ( 面 积 均 > 
70.00% ) ;强度 盐 渍 土 次 之 ,主要 分 布 在 腾 格 里 沙漠 
边缘 (人 研究 区 东北 - 北 - 南 沿线 边缘 ) 区 域 ,面积 随 土 
层 加 深 而 减 小 。 
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Spatial differentiation characteristics of soil salinity in Minqin Basin, 


downstream of Shiyang River, China 


LIU Xin, HAO Yuanyuan, HUA Limin 
(College of Patacultural Science, Gansu Agricultural University, Key Laboratory of Grassland Ecosystem, Ministry of 
Education; Engineering and Technology Research Center for Alpine Rodent Pest Control, National Forestry and 


Grassland Administration, Lanzhou 730070, Gansu, China) 


Abstract: Clarifying the spatial variability of soil salinity is the aim of soil salinization research and saline-alkali 
land improvements. In this study, the Mingin Basin in the lower reaches of the Shiyang River in the arid inland 
river basin was selected as the study area. By analyzing the statistical characteristics, spatial autocorrelation, 
variability, and distribution patterns of soil salinity in situ data of the 0-30 cm soil layer, the changes in the 
spatial heterogeneity of soil salinity could be explored to provide a theoretical basis for soil salinity in the study 
area. The measured soil salinity content and its variation range both decreased as the soil layer increased, and 
both were in a skewed normal distribution. Moran’ s Z index showed a “~>” type fluctuation trend, and the range 
of values increased with the soil layer, and the positive (negative) spatial autocorrelation decreased as the distance 
increased (negative spatial autocorrelation of 0- 10 cm was opposite). The exponential model was the most 
suitable for use with the variogram. The R’ increased with soil depth (0.62—0.69—0.81), and the RSS was 
0.0002. Variation rules of Co/(Co+C) and Ao were consistent, and the spatial heterogeneity ranged from weak to 
strong to moderate with soil depth. Accuracy of the interpolation was 10-20 cm > 20-30 cm > 0-10 cm. The soil 
salinity of the three layers had spatial differentiation characteristics, and the degree of salinization increased 
gradually from southwest to northeast in the horizontal direction (the same soil layer). In the vertical direction 
(different soil layers), the distribution patterns changed from simple to complex, and gradually tended to patches 
(10-20 cm) and spots (20-30 cm) from the aggregation state (0-10 cm). Nonsalinized soil (<1.50%) and four 
different levels of saline soil were distributed in the three soil layers, and the area of saline soil was moderate 
(>70.00%) > strength (approximately 20.00%) > mild (approximately 9.00%) > saline soil (<0.50 %). 


Keywords: soil salinity; spatial variation; geostatistics; Minqin Basin; downstrearm of Shiyang River 


